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Perfect-fluid  theory  is  applied  to  the  description  of  steady  motions  that 
can  be  generated  as  the  outflow  into  a  horizontal  channel  from  a  large  reservoir 
of  incompressible  heavy  fluid  whose  density  is  an  arbitrary  decreasing  function 
of  height.  A  particular  aim  is  to  pinpoint  the  significance  of  an  already  known 
class  of  flows,  called  self  similar,  which  satisfy  the  approximate  (shallow-water) 
equations  applicable  when  the  horizontal  scale  of  the  motion  greatly  exceeds  its 
vertical  scale,  but  which  have  not  until  now  been  shown  to  match  the  downstream 
conditions  that  primarily  determine  the  motion  in  practice. 

New  variational  principles  are  introduced  characterizing  the  class  of  self¬ 
similar  flows:  in  §2  there  is  a  characterization  in  terms  of  flow  force  among 
parallel  flows  realized  asymptotically  in  a  uniform  channel,  in  §3  among  a  wider 
range  of  possibilities  including  periodic  flows,  and  in  56  among  supercritical 
flows  realized  in  a  convergent-divergent  channel.  Aspects  of  general  flows  in 
channels  of  gradually  varying  breadth  are  treated  in  5 §4  and  5,  including  the 
remarkable  fact,  proven  in  §5,  that  every  steady  flow  outside  but  close  to  the 
self-similar  class  must  somewhere  undergo  a  local  crisis  unaccountable  by  the 
shallow-water  approximation.  Practical  interpretations  afforded  by  the  theoreti¬ 
cal  results  are  noted  in  §7.  - -  ~  S' 
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STEADY  FLOWS  DRAWN  FROM  A  STABLY  STRATIFIED  RESERVOIR 

.  .  i* 

T.  Brooke  Benjamin 

1.  INTRODUCTION 

It  has  been  known  for  some  time  that  the  equations 
describing  the  steady,  gradually  converging  flow  of  a  stably 
stratified  perfect  fluid  from  a  large  reservoir  into  a 
horizontal  channel  have  a  simple  solution,  according  to  which 
the  stream-surfaces  duplicate  those  in  a  corresponding  open- 
channel  flow  of  a  homogeneous  fluid.  Wood  (1968)  appears  to 
have  been  the  first  to  put  this  fact  on  record,  although 
evidently  it  was  then  already  known  to  others.  For  the  special 
flows  in  question  he  introduced  the  term  self  similar,  emphasizing 
their  property  that  the  heights  of  stream-surfaces  are  everywhere 
in  the  same  ratio,  and  this  useful  term  will  be  readopted  here. 

These  flows  were  also  demonstrated  by  Yih  (1969),  whose 
theoretical  description  went  a  little  further  than  Woods'  and 
included  an  appraisal  of  the  shallow-water  approximations  that 
underlie  the  mathematical  model.  Yih's  account  was  properly 
cautious,  moreover,  about  the  physical  significance  attributable 
to  the  special  class  of  flows.  Other  steady  flows  were 
recognized  to  be  possible  from  any  given  reservoir  that  either 
has  multiple  discrete  layers  or  is  continuously  stratified, 
and  it  was  noted  that  without  evidence  of  the  special  flows 
satisfying  terminal  conditions  that  can  be  imposed  downstream, 
there  is  no  a  priori  reason  for  their  ever  being  realizable. 

^Mathematical  Institute,  24/29  St.  Giles,  Oxford  0X1  3LB,  England 
Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 . 


Various  outstanding  questions  concerning  the  class  of 
self-similar  flows  will  be  answered  here,  most  notably  the 
question  of  their  realizability  in  practically  significant 
situations.  In  the  case  of  drainage  from  a  layer  in  the 
reservoir  bounded  by  a  discontinuity  in  density,  a  self-similar 
flow  will  be  shown  to  arise  as  the  extreme  state  before  hitherto 
stagnant  fluid  is  drawn  into  motion;  and  in  the  case  of  selective 
withdrawal  from  a  continuously  stratified  reservoir,  self¬ 
similar  flows  will  be  shown  always  to  have  priority.  The 
theoretical  models  studied  are  idealized,  talcing  no  account  of 
effects  due  to  viscosity  or  diffusion,  but  the  conclusions 
therefrom  should  have  bearing  on  certain  natural  flows  on  a 
large  scale  (cf.  Wood  1968,  $i) . 

Items  of  theory  are  developed  in  $$2-6,  and  finally  the 
powerful  interpretations  that  they  provide  are  assembled  in 
$7.  A  simple  but  very  informative  variational  principle  in 
terms  of  flow  force  (horizontal  pressure  force  plus  momentum 
flux)  is  established  in  $2  for  horizontal  flows  drawn  from  a 
reservoir  into  a  Straight  channel,  and  the  principle  is  extended 
in  $3  to  include  wavy  flows.  In  $4  the  theory  of  steady  flows 
in  a  channel  of  gradually  varying  breadth  is  reviewed,  the 
self-similar  solution  of  the  approximate  hydrodynamic  problem 
is  noted,  and  some  difficult  questions  posed  by  flows  other 
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than  the  self-similar  one  are  defined.  In  $5  it  is  proved 
that  no  steady  flow  neighbouring  on  the  self-similar  flow 
exists  corresponding  to  a  smooth  solution  of  the  shallow-water 
equations,  and  in  §6  a  second,  much  less  transparent  variational 

principle  is  given  relating  to  supercritical  flows  in  the 

* 

divergent  part  of  a  channel  with  a  throat. 

The  concluding  discussion  in  S7  relies  on  arguments 

that  generalize  the  following  elementary  idea  from  open-channel 

hydraulics.  Suppose  a  stream  of  water  having  asymptotic  depth 

h  and  velocity  u  is  drawn  from  a  large  reservoir  into  a  horizontal 

open  channel  of  rectangular  cross-section.  The  flow  force  of 

2  2 

t:he  stream  is  given  by  S  =  p(u  h  +  %  gh  )  ,  where  p  is  the  (uniform) 

2 

density?  and  the  Bernoulli  law  gives  u  »  2g(H  -  h),  where  H 

is  the  height  of  the  free  surface  in  the  reservoir  above  the 

3  2 

bottom  of  the  channel.  Thus  one  has  S  ■  pg(2Hh  -  -  h  ),  which 

2  2  2 

achieves  a  unique  maximum  value  S  «  —  pgH  when  h  ■  —  H  and 

m3  3 

2 

consequently  u  *»  gh.  If  the  reservoir  is  connected  through 
a  contraction,  the  flow  can,  of  course,  be  controlled  from 
the  downstream  end  of  the  channel,  say  by  lowering  a  weir  or 
by  operating  a  pump  which  takes  up  the  water.  In  every  case 
an  adjustment  increasing  the  steady  outflow  can  be  reckoned 
to  induce  an  increase  in  the  flow  force  of  the  oncoming  stream, 
which  has  to  match  the  greater  rate  of  extraction  of  momentum 
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from  the  system;  but  the  preceding  simple  result  shows  that 
upon  S  being  raised  to  the  value  S^,  no  further  increase  can 
in  any  way  be  induced.  It  is  impossible  to  realize  a  stream 
with  S  >  S  by  outflow  from  the  given  reservoir.  This 

IQ 

interpretation  of  the  critical  condition  at  which  downstream 
control  is  lost  may  be  appreciated  to  provide  better  physical 
insights,  and  to  admit  wider  generalization,  than  the  more 
usual  one  whifth  focuses  on  the  maximum  of  flow  rate  uh  for  a 
given  H.  An  immediate  advantage  of  it,  generalized  in  §3,  is 
that  the  maximum  principle  for  S  easily  extends  to  wavy  flows, 
every  one  of  which  is  known  to  realize  a  smaller  flow-force 
value  than  a  uniform  flow  with  the  same  H  (Benjamin  &  Lighthill 
1954). 

2.  STEADY  FLOWS  INTO  A  UNIFORM  CHANNEL 

As  indicated  in  figure  1,  a  stably  stratified  fluid 
lying  on  a  rigid  horizontal  plane  is  considered  to  flow  from 
an  infinitely  wide  reservoir  into  a  straight  channel,  whose 
cross-section  following  a  smooth  entry  region  is  rectangular  and 
uniform.  The  motion  is  assumed  to  be  steady,  and  the  fluid 
to  be  inviscid,  incompressible  and  non-dif fusive,  so  that  its 
density  has  a  constant  value  on  each  stream-surface,  varying 
only  among  these  surfaces.  Thus,  if  z  denotes  their  original 


Figure  1.  Illustration  of  flow  from  a  stratified  reservoir  into 
straight  channel. 


heights  above  the  bottom  in  the  reservoir  where  the  fluid  is 

at  rest,  the  density  is  representable  everywhere  by  o  “  p(z). 

All  pressures  will  be  stressed  relative  to  the  pressure 

at  the  level  z  ■  H  in  the  reservoir,  below  which  level  the 

fluid  is  drawn  into  the  channel  and  above  which  none  flows  when 

the  steady  motion  is  established.  According  to  the  hydrostatic 

law,  the  pressure  at  lower  levels  in  the  reservoir  is  given  by 

H 

p  -  g  J  pU)  d?. 


and  hence  the  total  head  of  the  fluid  above  the  bottom  by 


R(z)  -  gzp  +  p  -  g{zp  +  [  p(§)  d§) 

Jz 

H 

-  g{psH  -  J  Sp'U)  <H), 


(2.1) 


where  pg  =  p(H-).  We  also  need  the  notation  p+  =*  p(H+), 
supposing  for  now  that  p  <  p  (see  figure  1) ,  but  the  case 
p.  »  p  will  be  covered.  In  (2.1),  p'  denotes  dp/dz,  being 

»  D 

a  non-positive  function  by  the  assumption  of  stability,  and 
henceforth  accents  will  be  used  only  to  denote  derivatives  of 
z-dependent  functions. 

If  the  flow  approached  asymptotically  along  the  channel 


is  horizontal  (i.e.  free  from  waves),  the  height  y  of  the 
stream-surfaces  in  it  is  a  function  of  z  alone,  and  we  write 
y(H-)  *  h.  Since  the  fluid  is  at  rest  at  heights  above  h. 


the  layer  between  heights  h  and  H  must  be  filled  with  stagnant 


fluid  of  density  p+.  Hence  the  pressure  in  the  horizontally 
moving  fluid  is  given,  according  to  the  hydrostatic  law,  by 

h  H 

p  a  g{p  (H  -  h)  +  f  p  dy>  =  g{ p  (H  -  h)  +  /  p(5)y’(5)d£} 

y  z 

H 

-  g{p+(H  -  h)  +  Psh  -  py  -  /  yp'  d£}.  (2.2) 


If  q  is  the  magnitude  of  the  velocity  in  the  fluid,  the 
Bernoulli  law  tells  us  that 

2 

R  »  gyp  +  p  +  *Sp  q  (2.3) 

is  constant  on  any  stream-surface,  specified  by  the  value  of 

2  2 

z;  and  in  the  asymptotic  flow  q  *  u  ,  the  square  of  the 
single,  horizontal  component  u(z)  of  velocity.  Hence,  using 
the  expressions  (2.1)  for  R(z)  and  (2.2)  for  p,  we  obtain 

2  A  H 

pu  «  2g{ p (H  -  h )-/(£-  y)  p 1  d£}  ,  (2.4) 

z 

in  which  p  a  .  p+. 

The  key  to  subsequent  physical  interpretations  is  to 

consider  the  flow  force  S  defined  as  the  sum  of  horizontal 

pressure  force  and  momentum  flux,  per  unit  span,  in  the  whole 

layer  of  fluid  affected  by  the  steady  flow  (cf.  Benjamin  1966,  §2) 

2 

Thus  S  is  the  integral  of  p  +  pu  with  respect  to  height,  from 
the  bottom  to  the  height  H.  The  contribution  to  S  from  the 
pressure  in  the  stationary  fluid  above  the  uppermost  stream- 
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ft, 


w 


2 

surface  y  ■  h  is  plainly  %gp+(H-h)  ,  and  hence 

9  h  2 

S  ■  %gp  (H  -  h)  +  /  (p  +  pu  )dy 

+  o 

JJ 

-  -  h)2  +  /  (P  +  pu2)y '  dz.  (2.5) 

+  o 

For  substitution  in  this  integral,  p  can  be  expressed  by  (2.2) 

2  2  2 

and  the  momentum  flux  density  pu  by  (2.4)  with  q  =»  u  . 

We  thus  obtain  at  once 


S/g  a  *5p  (H  -  h)  2  +  p  (H  -  h)h  +  P  h2  +  2p  (H  -  h)h 

*  T  S 

H  H 

-  J  (pY  -  J  [2&  -  y(§)]  p*  (SJdgJy'dz, 


o  z 

and  after  an  integration  by  parts 
.2  .  a _  3  .2 


H 


S/g  -  >5p+H2  +  p(2Hh  -  |  h2)  +  j  (-P')  (2zy  -  |  y2)dz.  (2.6) 

o 

With  regard  to  the  terms  of  (2.6),  note  the  identity 


„  3  2  2  2  1  _  ,2 

2zy  -  -  y  »  -  z  -  -  (2z  -  3y)  , 


including  its  instance  with  z  *  H,  y  =  h.  Note  also,  from  the 
assumption  of  stable  stratification  in  the  reservoir,  that 
p  >0  and  -  p'  is  a  non-negative  function  on  [0,H) .  It 
therefore  follows  from  (2.6)  that  S  is  an  absolute  maximum 
when 


y 


(2.7) 


which  case  is  the  self-similar  flow  that  is  critical  in  a 
sense  to  be  recalled  presently. 
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The  maximum  value  of  S  is 


r ,  „2  2  A  2  2c.  2  .  , 

Sm  *  ««P+H  +  J  PH  +  J  j  (-  p')z  dz} 

o 

H 

-  SQ  +  i  g{J5pgH  +  j*  (-  p‘)z  dz}. 


(2.8) 


2  |»  o 

wljere  S  a  g{%p  H  +  %  |  (-  p’)z  dz} 

O  S  J 

o 

H 

*  g  J  pz  dz 


is  the  value  of  S  for  the  fluid  in  a  state  of  rest.  This 
result  generalizes  the  well-known  principle  recalled  in  SI 
concerning  open-channel  flows  of  a  homogeneous  fluid,  a  case 
that  is  recovered  by  putting  p+  a  0  and  p*  »  0  in  (0,H). 

For  the  self-similar  flow  described  by  (2.7),  we  have 
from  (2.4)  that 


pu2  »  j  g{pH  -  J  o'z  dz}  =  |  g{R(z)  -  R(H+)} 


h  ~ 


g(ph  -  J  Y  dy}. 


(2.9) 


(2.9’) 


where  p  a  ”p(y)  53  ?K(jz)  is  written  in  place  of  p(z).  From 
the  definition  u  a  d^/dy,  *(0)  =  0,  (2. 9')  can  be  used  to  find 
the  stream-function  ${y)  for  this  flow. 
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A  convenient  example 


The  example  now  introduced  will  be  used  later  to 
illustrate  other  aspects  of  the  problem.  Let 


H  ■  1,  p(z)  =*  { 1  -  pz)  (0  <  P  <  1)  , 

ps  *  po(1  ”  ^  =  p+  (i*e*  P  =  0)> 


(2.10) 


In  this  case,  (2.9)  gives 


1  o 

The  velocity  is  thus  a  maximum  at  y  =  jB  +  0(B  )  (i.e.  just 

above  the  bottom  of  the  channel  if  3  is  small)  and  falls 
continuously  to  zero  as  y  -»  j. 


3.  GENERALIZATION  OF  THE  VARIATIONAL  PRINCIPLE 

It  has  been  shown  that  among  all  horizontal  steady  flows 
realizable  by  withdrawing  fluid  from  a  bottom  layer  of  a 
reservoir  with  any  given  density  stratification,  the  self¬ 
similar  flow  described  by  (2.7)  uniquely  achieves  the  maximum 
possible  flow  force  S.  This  principle  will  now  be  extended 
to  two-dimensional  wavy  flows,  which  plainly  may  be  realized 
in  certain  circumstances,  as  when  the  flow  passes  into  a  uniform 


stretch  of  channel  downstream  of  an  obstacle. 


4rfSw>--- 


Por  any  two-dimensional  flow  of  a  stratified  but 
incompressible  fluid,  in  which  the  velocity  components  are 
u  and  v  respective  to  x  and  y,  there  is  a  stream- function  t 
by  virtue  of  the  fact  that  div(u,v)  ■  0.  Moreover,  when  the 
flow  is  steady,  the  streamlines  $  ■  const,  coincide  with 
the  lines  p  *  const,  and  so  also,  in  the  present  description, 
with  the  lines  z  *  const.  We  may  accordingly  take  y  *  y(x,z) 
as  the  dependent  variable  and  use  the  fact  that  f,  p  and  R 
are  functions  of  z  alone.  The  velocity  components  are  then 
given  by 

U  -  ty  -  ♦'<*>  y  ' 

yx 

v  ■  “  “  ♦'  (z)7i  * 

x  yz 

and  it  is  helpful  to  write 

F(z)  -  p(z)  [*’  (z)]2. 

In  general  this  function  is  not  prescribable  a  priori,  but 

for  all  self-similar  flows  (see  §4)  it  is  proportional  to 

the  function  of  z  expressed  on  the  right-hand  side  of  (2.9). 

On  the  assumption  that  y  >0,  the  partial  differential 

z 

equation  satisfied  by  y(x,z)  may  be  found  easily  enough  from 
the  Euler  equations  of  steady  motion,  or  even  more  readily 
by  transforming  the  equation  introduced  by  Long  (1953)  for 

le 

♦(x,y).  as  simplified  by  Yih  (1960)  as  an  equation  for  p  $. 


-u- 


It  is 


-  F(zHyx/yzlx  +  %[?(«)  d  +  yx)/yz32 


+  gp* (z)y  -  R* (z)  «  0.  (3.1) 

Except  that  the  coefficient  function  F  is  undetermined,  this 
equation  resembles  one  that  has  been  used  previously  for 
somewhat  simpler  problems  of  internal  waves  assumed  to  arise 
from  a  base  flow  with  uniform  horizontal  velocity  (cf.  Benjamin 
1967,  eqn.  (3.3),  also  Turner  1960  where  an  exacting  mathematical 
treatment  is  given).  The  boundary  conditions  to  be  satisfied 
by  the  solution  y  of  (3.1)  are  the  kinematical  condition  at 
the  bottom  of  the  moving  fluid, 

y(x,0)  =  0  V  X,  (3.2) 


and  the  dynamical  condition  ensuring  that  at  the  top  z  =*  H-, 

with  y(x,H-)  =  h(x),  say,  the  pressure  is  p+(H  -  h) .  When 

this  pressure  is  alternatively  expressed  by  (2.3)  with 

R(H-)  3  gp  H,  the  condition  is  seen  to  take  the  form 
s 


1  +  h 


*5F(H) 


gp(H  -  h)  V  x. 


(3.3) 


In  the  case  that  p  =  0  and  consequently  F  i  0  as  z  f  H  (see 

example  in  §5),  the  upper  boundary  condition  becomes  simply 

that  y  remain  bounded  in  the  limit, 
z 
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For  present  purposes  it  is  sufficient  to  know  that  any 
steady  wavy  flow  must  he  represented  by  a  solution  of  (3.1) 


•satisfying  these  boundary  conditions,  and  our  inability  to 

prescribe  the  function  7  is  in  fact  immaterial.  The  flow 

force  S  is  again  given  by  (2.5).  from  which  p  can  be  eliminated 

2  2  2 

by  means  of  (2.3)  with  q  •  u  +  v  . 

The  result  is 

2  rBr  '  V 

S  -  *SP+( H  -  hr  +  J  {%P(z)(  —  x)  +  Rys  -  gpyyz)  dz, 

o 

which  an  integration  by  parts  with  use  of  (2.1)  and  (2.3) 
reduces  to 

S  «  *sgp+  H2  +  gp(Hh  -  %h2) 

*  2 

H  1  -  y 

+  f  {*SF(z)( — - — -  gp*  (2)  (zy  -  Hy  )}  dz.  (3.4) 

By  differentiation  of  this  expression  and  by  appeal  to  the 
fact  that  y  satisfies  (3.1)  and  the  boundary  conditions  (3.2) 
and  (3.3).  it  is  easy  to  confirm  that  S  is  independent  of  x, 
as  is  plainly  required  by  momentum  conservation.  Equation  (3.1) 
and  the  boundary  conditions  can  be  used  further  to  reduce  the 
terms  in  (3.4)  dependent  on  F(z).  Integrations  by  parts  lead 
directly  to 
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mgjjgjg 


mdHHsinti  ii 


*• —  • ••  • 


(3.5) 


S/g  -  Jjp+  H2  +  p(2Hh  -  |  h2) 


+  J  (— p * )  (2zy  -  ly2)dz  -  J  F(z)|—  +  y(~)x}  dz 


which  recovers  (2.6)  in  the  case  of  x-independent  flows. 
Furthermore,  since  S  is  independent  of  x  and  since  the  average 
of  the  second  integral  in  (3.5)  is  zero  over  a  wavelength 
of  a  flow  periodic  in  x  (or,  more  generally,  is  zero  between 
any  two  stations  where  y  is  the  same  or  where  y^  »  0  V  z  6  (0,H)), 
we  can  conclude  from  (3.5)  that  for  a  wavy  flow 


S  m average  of  SQ(y),  (3.6) 

where  S  (y)  denotes  the  expression  (2.6)  with  x-dependence 
o 

of  y  now  allowed. 

In  the  light  of  the  discussion  following  (2.6),  it  follows 

immediately  from  (3.6)  that  for  every  wavy  flow  S  is  less 

than  the  value  S  realized  uniquely  by  the  special  flow  that 

m 

(2.7)  describes.  The  promised  generalization  of  the  variational 
principle  is  thus  established. 

Note  also  that  whatever  the  average  y(z)  of  y(x,z) 
with  respect  to  x  in  a  wavy  flow,  (3.6)  implies  that  S  <  SQ(y). 
This  conclusion  accords  with  a  property  well  known  from  studies 
of  other  formulations  of  internal-wave  problems,  namely  that 
when  periodic  waves  can  be  superposed  without  energy  loss  on 
an  originally  horizontal  flow,  the  process  entails  a  reduction 
in  flow  force  (cf.  Benjamin  1966,  §2) . 
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4.  GRADUALLY  VARYING  FLOWS 


Let  ua  extend  the  ideas  of  $2  to  the  case  of  stead/ 

flow  from  a  reservoir  into  a  channel  with  planform  as 

illustrated  in  figure  2.  The  breadth  b  of  the  channel  varies 

continuously  with  horizontal  distance  x,  having  a  minimum 

value  b  as  indicated  in  the  figure,  and  its  variation  is 
c 

assumed  to  be  so  gradual  that  the  velocity  u  of  the  fluid  in 
the  x-direction  is  the  only  component  significantly  entering 
the  Bernoulli  law  (2.3).  Moreover,  u  is  taken  to  be  uniform 
across  the  span  of  the  flow,  although  of  course  its  value 
depends  on  height  above  the  bottom.  The  theory  to  be  developed 
is  thus  a  counterpart  of  the  shallow-water  approximation  for 
gradually  varying  open-channel  flows  of  a  homogeneous  fluid. 

Conservation  of  mass  in  the  incompressible  fluid  can 
be  expressed  by  considering  an  elementary  stratum  of  vertical 
thickness  6z  in  the  reservoir,  from  which  the  fluid  flows 
into  a  sheet  whose  local  thickness  in  the  channel  is  fiy. 

If  u  is  locally  the  mean  horizontal  velocity  in  the  sheet,  then 

ubjy  ■  &Q 

is  the  element  of  volume  flux  which  must  be  independent  of 
x  when  the  flow  is  steady.  Dividing  by  ftz  and  taking  the 
limit  as  ftz  +  0,  we  have 
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-i*a 


uby*  -  Q', 


(4.1) 


where  the  s-derivative  Q*  is  evidently  a  function  of  z  alone. 

2 

The  expseeeion  (2.4)  for  pu  is  applicable  ae  an  approximation 
everywhere  according  to  the  present  assumptions  (cf .  Yih  1969) , 
and  the  elimination  of  u  between  it  and  (4.1)  gives 

2  a  « 

2Xy'  { p(H  -  h)  -  f  (5  -  y)p'(5)d5}  -  f(z>,  (4.2) 

"z 

in  which  X  ■  b2/b2,  h  ■  y(H)  and 

c 

f(s)  -  p(z) [Q • (z) ] 2/gb2  . 

This  equation  for  y(z)  on  (0,H)  is  complemented  by  the 
boundary  condition  y(O)  ■  0,  and  the  dependence  of  the  flow 
on  its  position  along  the  channel  enters  through  the  parameter 
X  which  varies  from  1  to  •. 

The  problem  may  alternatively  be  expressed  as  a  second- 

order  differential  equation  with  a  pair  of  boundary  conditions. 

2 

Dividing  (4.2)  by  y'  and  differentiating  with  respect  to  z, 
we  obtain 

(f/y’2)'  +  2X(y  -  z)p*  -  0,  (4.3) 

and  the  boundary  conditions  are 

y<0,  x)  -  0,  2Xp{H  -  y(H,\)}[y'  (H,X)2  -  f(H-).  (4.4) 

This  form  of  the  problem  corresponds,  of  course,  to  the 
x-independent  version  of  (3.1)  -  (3.3)  with  F(z)  ■  gf(z)/x. 

As  noted  earlier,  in  the  case  that  p  »  0  and  consequently 
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f (H— )  m  0,  the  second  condition  in  <4. 4)  is  replaced  by  a 
condition  of  regularity  as  z  •»  H,  i.e.  ly'(H-,\)|  <  •. 

The  non-linear  parametrized  problem  (4.3)  and  (4.4) 
is  not  amenable  to  comprehensive  treatment  that  allows  arbitrary 
specification  of  the  coefficient  function  f.  A  somewhat 
perplexing  situation  is  presented,  moreover,  in  that  f  may 
not  be  a  prescribable  feature  of  the  physical  problem.  It 
will  generally  depend  on  conditions  imposed  at  the  downstream 
termination  of  the  flow,  and  there  appears  to  be  no  simple 
argument  delimiting  the  complete  class  of  functions  f  that 
are  physically  relevant.  However,  by  analogy  with  the 
corresponding  problem  for  open-channel  flows  of  a  homogeneous 
fluid,  we  may  reasonably  look  to  the  possibility  of  solutions 
with  the  following  properties: 

(i)  y(z,\)  varies  continuously  with  \  >  1  (hence 

continuously  with  distance  x  along  the  channel) . 

(ii)  The  solution  in  1  <  \  <  »  has  two  branches  which 
are  confluent  at  X.  *  1  (i.e.  the  flow  can  be  different  at 
channel  sections  with  the  same  breadth  upstream  and  downstream 
of  the  minimum  section). 

(iii)  On  one  branch,  y(z,X)  z  as  \  -»  •  (i.e.  the  flow 
connects  smoothly  with  the  reservoir) . 


Self-similar  flows 


As  was  appreciated  in  §1*  these  special  flows  have 
previously  been  noticed  to  be  solutions  of  the  above  problem. 

If  we  take 

y  -  *(\)z,  (4.5) 

which  satisfies  the  boundary  condition  on  the  channel  bottom, 
then  equation  (4.1)  becomes 
3  R 

2X(h  -  h  )  { oH  -  [  So'(S)  <3g)  -  f(z). 

Jz 

Hence  (4.5)  is  a  solution  with  the  required  properties  if 

X(*2  -  H3)  ■  ^  (4.6) 

and 

f<z)  -  ~(oH  -  J  So'(g)  dg}  -  (4.7) 

2 

2 

The  cubic  (4.6)  for  *  has  a  double  root  h  «■  —  for 

c  3 

X  ■  1;  and  it  has  two  distinct  positive  roots  for  X  >  1,  one 

of  which  tends  to  1  and  the  other  to  zero  as  \  «♦  ».  Thus 

the  properties  (i),  (ii)  and  (iii)  are  all  provided.  In  fact, 

h  «  x(X)H  is  precisely  the  local  depth  of  a  ’choked'  flow 

of  homogeneous  heavy  fluid  along  an  open  channel  of  the  same 

planform  \  »  x(x)  (i.e.  the  flow  through  a  Venturi  flume, 

2 

which  is  subcritical  (u  <  gh)  upstream  and  supercritical 
downstream  from  the  minimum  sectionj. 
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It  remains,  of  course,  to  demonstrate  the  significance 


of  this  simple  result.  Since  other  flows  are  generally 
possible  from  stratified  reservoirs,  there  is  a  need  to 
explain  how  and  why  self-similar  flows  might  be  generated. 


Regular ity  of  solutions 


Let  us  examifte  conditions  for  a  solution  y[z,X(x)] 
of  (4.3)  and  (4.4)  to  vary  smoothly  with  distance  x  along 
the  channel.  Assuming  \(x)  to  be  continuously  differentiable, 
differentiating  (4.3)  with  respect  to  x  and  writing 


qp(x,\) 


2L  m  SSL 
ax  dx  3X  ' 


-we  obtain 


~  Xp'cp  -  (y  -  z)p'  •  (4*8) 

Similarly,  the  boundary  conditions  for  <p  are  seen  from  (4.4) 
to  be 


<p(0)  =  0,  y'(H)tp(H)  -  2(H-h)?'(H)  =  .  (4.9) 


(Here  the  dependence  of  ®  and  y  on  X  is  left  implicit.)  For 
a  smooth  solution  of  the  parametric  problem  (4.3)  and  (4,4), 
the  linear  problem  (4.8)  and  (4.9)  for  must  also  have  a 
solution  for  each  relevant  value  of  the  parameter  X. 
Accordingly,  although  the  issue  is  trivial  in  the  case  of 
self-similar  flows,  something  about  the  general  case  may  be 


learned  by  means  of  the  Fredholm  alternative  principle. 

Inquiries  on  this  basis  have  so  far  made  limited 
progress.  Some  comparatively  easy  conclusions  will  be  noted 
below,  but  the  main  proposition  in  view  can  only  be  stated  as 
a  conjecture,  made  plausible  by  the  outcome  of  the  linearized 
perturbation  theory  to  be  presented  in  §5.  It  is  that  except 
for  the  self-similar  flow  described  by  (4.5)  -  (4.7),  the 
system  (4.3)  and  (4.4)  has  no  continuous  solution  extending 
to  the  minimum  section  where  X  «  1,  dx/dx  «  0,  and  being 
such  that  qj  fc  0  there. 

At  the  minimum  section,  a  possibility  according  to 
(4-8)  and  (4.9)  is  that  cp  ■  0,  in  which  case  the  solution 
y(z,x)  on  the  downstream  side  returns  along  the  same  branch 
as  it  approaches  X  «  1  on  the  upstream  side.  For  reasons  that 
will  be  noted  in  §7,  this  possibility  has  comparatively  little 
interest,  and  it  will  not  be  considered  further  here.  Thus, 
at  X  =  1 ,  cp  is  required  to  be  a  non-trivial  solution  of  the 
homogeneous  version  of  the  linear  boundary-value  problem  (4.8) 
and  (4.9). 

The  meaning  of  this  requirement  is  made  clearer  by 
transforming  the  left-hand  sides  of  (4.8)  and  (4.9)  so  that 
the  local  height  y  of  the  stream  surfaces  is  the  independent 
variable.  Thus,  writing  cp  =  ^(y)  and  using  (4.2)  together 
with  (2.4),  one  derives 


(pu  ny)y  -  gpn  =  0, 

2 

n(0)  =  0,  pu  ny  =  gpn  at  y  =  h. 


The  existence  of  a  non-trivial  solution  of  (4.10)  is 
recognizable  as  being  just  the  condition  for  an  infinitesimal 
wave  of  extreme  length  to  be  superposable  without  energy  loss 
on  a  horizontal  flow  that  has  the  given  velocity  u  =  u(y) 
and  density  distribution  p  =  p (y)  [cf.  Benjamin  1966,  §3.3] 
So  the  flow  at  X  =  1  is  critical  in  the  usual  sense  of  the 
term.  This  result  is  to  be  expected,  of  course,  because 
under  the  present  shallow-water  assumptions  the  condition 
that  the  flow  have  a  smooth  non-zero  variation  through  the 
throat  of  the  channel  is  evidently  equivalent  to  a  horizontal 
flow  at  A  =  1  admitting  a  long-wave  perturbation. 

It  may  straightforwardly  be  shown  that,  in  all  cases, 
the  accelerating  flow  in  the  divergent  part  of  the  channel 
is  supercritical  in  the  sense  that  a  free  infinitesimal  long 
wave  propagating  against  it  would  be  swept  downstream.  This 
property  corresponds  to  \  <  where  Yj  is  the  lowest 

eigenvalue  of  the  homogeneous  Sturm-Liouville  problem  related 
to  (4.8)  and  (4.9)  [i.e.  with  y  in  place  of  X  and  dX/dx 

The  existence  of  <9  everywhere  downstream  of  the  throat  is 
hence  ensured  according  to  the  Fredholm  principle. 


5.  FLOWS  CLOSE  TO  THE  SELF-SIMILAR 


Again  with  regard  to  a  channel  of  gradually  varying 
breadth,  we  consider  steady  flows  that  are  small  perturbations 
from  the  self-similar  flow  defined  by  (4,5)  -  (4.7),  The 
solution  is  now  expressed  in  the  form 

y  »  kz  +  c£(z,k) ,  (5.1) 

where  the  parameter  k  varies  from  0+  to  1,  being  related 
to  X  by  (4.6),  and  t  is  an  infinitesimal  number.  For 
simplicity  of  illustration  the  discussion  refers  particularly 
to  the  example  specified  in  (2.10),  but  the  means  will  be 
indicated  whereby  the  conclusions  can  readily  be  extended 
to  all  other  examples. 

In  the  chosen  example,  equation  (4.3)  becomes 

(o/y'2)'  -  2\(y  -  z)  *=  0  (5.2) 

on  0  <  z  <  1,  with  a  «  f(z)/pQ0,  o(l)  **  0;  and  the  boundary 
conditions  are 

y(0,n)  =»  0,  ly'  (1,k)I  <  »  •  (5.3) 

In  keeping  with  (5.1),  the  coefficient  function  c(z)  is 
represented  as  an  infinitesimal  perturbation  from  its  form 
given  by  (4.5)  for  the  self-similar  flow,  thus 

o(z)  =  (1  -  z2)  +  cr(z).  (5.4) 
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The  perturbed  flow  is  assumed  to  originate  from  the  same 
bottom  layer  in  the  reservoir,  and  so  it  is  implied  that 
r  (1)  =  0. 

After  substitution  of  (5,1)  and  (5.4)  into  (5.2), 
linearization  in  e  gives 

{(1  -  z2) £' } '  +  =  ar*,  (5.5) 

in  which  |j  =»  (27/4)  Xx'*  *  x/(l-x) 

and  a  »  (27/8) x 

have,  like  x,  two  positive  values  for  each  X  >  1.  The 

2 

number  jj  equals  2  at  the  critical  section  (x  =  ~)  and 

increases  smoothly  with  X  on  the  upstream  side,  with  jj  -♦  » 

2 

-as  X  +  ®,  x  t  1.  On  the  downstream  side,  k  <  —  and  therefore 

U  <  2.  The  required  solution  £  of  (5.5)  must,  of  course, 

satisfy  the  boundary  conditions  (5.3). 

Now,  the  homogeneous  equation  corresponding  to  (5.5) 

is  Legendre's  equation,  which  has  a  non-trivial  solution 

bounded  on  [-1,1],  the  respective  Legendre  polynomial  P  (z) , 

m 

when  jjL  takes  the  succession  of  values  m(m+l)  (m  =  1,2,  ...). 

In  view  of  the  first  boundary  condition,  only  the  odd-order 

polynomials  which  vanish  at  z  =  0  are  relevant  here,  but 

2 

they  comprise  a  basis  in  L  (0,1).  Thus  it  is  merely  enough 
2 

that  r*  e  L  for  there  to  exist  a  representation 
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r'  ‘  S  %  *2n-l,z>‘  an  ‘  (4n-l)  j  r'  p2n-l(z)  dz<  (5<6) 
n«l  o 

2 

in  which  the  sequence  of  coefficients  {a^}  e  f  .  Accordingly, 
the  formal  solution  of  (5.5)  satisfying  the  boundary 
conditions  is 


C 


CL 

"n^l  H  "  2n(2n~1) 


a 

n 


2n-l 


(z)  , 


(5.7) 


which,  since  a  is  a  bounded  positive  number,  is  meaningful 

except  where  the  denominator  jj,  -  2n(2n  -  1)  vanishes  for  any 

coefficient  with  a^  0.  [Note  incidentally  that,  with  these 

2 

exceptions,  the  attribution  r*  e  L  implies  the  existence 

of  at  least  a  weak  solution  £  e  H^(0,1),  and  further  regularity 

of  £  follows  from  that  or  r.  For  example,  if  r  e  C1, 

2 

then  £  e  C  and  the  ordinary  differential  equation  (5.5) 
is  satisfied  pointwise.J 

The  first  conclusion  to  be  drawn  from  (5.7)  is  that  if 
the  perturbed  flow  is  to  exist  at  a  throat  in  the  channel, 
where  \i  -  2,  then  a^  =  0  in  the  expansion  (5,6).  On  the 
downstream  side,  we  have  that  |i  -  2n(2n-l)  <  0  for  all 
n  =  1,2,...,  and  thus  the  solution  (5.7)  remains  meaningful 
everywhere  in  the  supercritical  region.  On  the  upstream  side, 
however,  every  one  of  the  numbers  2n(2n-l)  >  12  is  crossed 
by  ji,  and  so  (5.7)  is  valid  everywhere  only  if  a^  =  0  for 
all  n.  The  completeness  of  the  basis  {P„  ,}  hence  establishes 

the  significant  conclusion  that  no  smooth  solution  of  the 
shallow^water  equations  exists  neighbouring  on  the  sel f-s ini lar 
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flow  from  the  given  reservoir.  Every  perturbed  steady  flow 
suffers  at  least  one  lbcal  crisis  where  the  shallow-water 
approximation  ceases  to  be  valid,  however  gradual  the  variation 
in  breadth  along  the  channel. 

This  conclusion  readily  extends  to  all  other  examples 

of  stable  density  stratification  in  the  reservoir.  For  each, 

in  place  of  (5.5),  an  equation  with  another  Sturm-Liouville 

operator  on  the  left-hand  side  will  be  posed,  and  if  p  >  0 

the  upper  boundary  condition  will  be  given  by  linearizing 

the  second  of  (4.4) .  By  the  Riesz-Fischer  theorem,  the  set 

of  eigensolutions  for  each  respective  Sturm-Liouville  problem 

2 

is  complete  as  a  basis  in  L  (0,H),  and  accordingly  the 
argument  proceeds  as  above. 

It  deserves  emphasis  that  a  failure  of  the  shallow- 
water  approximation  in  some  part  of  the  flow  does  not  necessarily 
invalidate  the  approximate  solution  (5.7)  elsewhere.  The 
nature  of  the  local  crises  that  all  flows  neighbouring  the 
self-similar  have  been  shown  to  suffer  will  not  be  explored 
in  any  detail  here,  although  it  is  an  interesting  matter  that 
should  be  worth  further  study.  The  phenomena  indicated  are 
presumably  continuous  processes  according  to  a  more  accurate 
perfect-fluid  model,  but  to  describe  them  one  needs  to  abandon 
the  hydrostatic  approximation  for  pressure  and  use  explicitly 
x-dependent  differential  equations  such  as  (3.1).  A  particular 
component  with  n  >  1  in  the  Sturm-Liouville  expansion  of 
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C(z,h)  will  be  better  modelled  in  the  vicinity  of  its  crisis 

point,  say  x  ■  0,  by  a  function  of  the  form  w^(x)5n(z),  where 

5n(x)  is  the  eigenfunction  in  question  (e.g.  P2n-1^  above) 

and  where,  to  a  first  approximation,  a  suitably  normalized 

w  (x)  will  satisfy  the  equation 
n 

k*"2(d2w  /dx2)  -  xw  »  1, 
n  n 

-2 

in  which  k  is  a  positive  parameter  proportional  to 

-2 

dy/dx  >  0  at  x  *  0.  The  specification  k  *  0  recovers 

the  unacceptable  local  singularity  given  by  the  shallow-water 

-2 

approximation.  But  with  k  >0  this  equation  has  a  solution 
that  is  bounded  on  -  «  <  x  <  ».  Expressible  in  terms  of 
Airy  functions,  the  needed  solution  is  -  tj  Gi(kx)  in  the 
notation  used  by  Abramowitz  &  Stegun  (1965,  p.448).  For  large 
kx  >  0,  the  solution  is  quickly  asymptotic  to  -  x  1;  and  for 
large  kx  <  0, 

w  ^  ir^-kx)  cos  [--(-kx) 3^2  +  7tt]  -  x  1  +  0[(-kx)  7/^4] . 
n  o  •* 

The  slowly  diminishing  but  increasingly  rapid  oscillations 

may  roughly  simulate  what  in  fact  happens  upstream  of  a  crisis 


6.  THE  EXTREMAL  PROPERTY  OP  SUPERCRITICAL  SELF-SIMILAR  FLOWS 

We  return  to  the  flow-force  principle  demonstrated  in 

S2,  now  considering  it  subject  to  the  constraint  that  the 

flows  competing  for  the  maximum  of  S  have  passed  through  the 

throat  into  the  divergent  part  of  the  channel  and  so  have 

become  supercritical.  Let  subscript  1  refer  to  any  particular 

2 

station  downstream  (i.e.  >1,  k.  <  — ) ,  and  let  S  denote 

x  x  3  In 

the  value  of  S  given  by  substitution  of  y  »  in  (2.6). 

Taking  the  expression  (5.1)  for  y  but  no  longer  assuming  e 

to  be  infinitesimal,  we  obtain  from  (2.6) 

H 

S±/g  -  Slm/g  +  e(2-3K1){pHC1(H)  +  J  (-o') zqdz} 

o  ~ 

-  f  €2{pq(H)  +  jVp'K*  dz}.  (6.1) 

O 

Here  the  coefficient  of  c  is  to  be  treated  as  the  first 
variation  S^/g  of  S^/g  among  flows  differing  from  the  self¬ 
similar  flow. 

On  the  other  hand,  when  (5.1)  is  substituted  and  (4.6) 

is  used,  (4.2)  leads  to 

H  H 

r  «  J  %f(z)  dz  =  (4/27) { pH2  -  /  p'z2  dz} 

o  o 

„  H 

+  c\k(2  -  3n) { pH? (H, p)  -  /  p'Cz  dz} 

z 

2  H  ,  A  H 

+  e  X  [  <1  -  n )/  V  (pH  -  /  p'(l)C  <3U 
o  z 

H 

-  2H{pC (H,h)  -  J  p*c  dz}]  +  OU  ), 


(6.2) 


which  by  the  definition  of  f(z)  must  be  the  same  for  all  k 

such  that  y(z,x)  exists.  Note  that  whereas  local  failures 

of  the  shallow-water  approximation  may  occur  upstream  of  the 

throat,  (6.2)  still  holds  everywhere  downstream.  Writing 

(6.2)  as  r  ■  r  +er+*5er+...,we  have  that  f*  «.0  at 
o 

2 

the  throat  where  h  **  j  and  therefore,  in  the  limit  e  ■»  0, 

.  2 
also  r  ®  0  everywhere  downstream  where  0  <  x  <  ”  .  Since 

*1*1^1  *  gr±  according  to  (6.1)  and  (6.2),  it  follows  that 

-  0.  (6.3) 

Thus  S  is  a  stationary  value  for  variations  in  the  class 
lm 

of  flows  that  are  smooth  at  and  downstream  of  the  throat. 

The  study  of  second  and  higher  variations  is  more 
complicated  and  will  merely  be  outlined  here.  One  proceeds 
by  using  (6.2),  in  the  form  of  the  identity 

er(x±)  =*  Jse2{r(|’)  -  rU^))  +  o(e3), 

to  reduce  the  e-term  in  (6.1).  The  results  of  the  linearized 

theory  summarized  in  §5  are  then  used  to  evaluate  S  .  It 

is  thus  shown  without  much  difficulty  that  0.  Further 

estimates  of  S  -  S  finally  confirm  that  is  a  maximum. 

1  lm  lm 
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7.  PHYSICAL  CONCLUSIONS 


Bven  within  the  context  of  perfect-fluid  theory,  the 
general  problem  of  selective  withdrawal  from  a  stratified 
reservoir  is  largely  intractable  because  of  the  freedom 
evidently  available  in  posing  the  downstream  conditions  that 
determine  the  flow.  While  falling  far  short  of  a  general 
solution,  the  preceding  results  nevertheless  illuminate 
various  apparently  central  aspects  of  the  problem,  and  on 
the  basis  of  them  a  number  of  plausible  interpretations  can 
be  made  as  follows  about  practical  possibilities. 

1.  First  take  the  case  of  stratified  fluid  drawn  steadily 
into  a  straight  channel,  as  illustrated  in  figure  1,  and 
suppose  that  the  flow  is  caused  by  the  extraction  of  fluid 
through  a  slot  near  the  bottom  at  the  end  of  the  channel.  If 
there  is  a  discontinuity  of  density  at  z  *  H  (i.e.  p  >  0) , 
it  is  intuitively  clear  that  the  fluid  above  the  interface 
will  not  be  drawn  into  the  slot  when  the  extraction  rate  is 
sufficiently  small,  and  so  a  definite  question  arises  about 
the  limiting  condition  beyond  which  the  fluid  originally  above 
z  «  H  in  the  reservoir  begins  to  be  extracted.  In  other  words, 
what  is  the  'drawdown  condition'  corresponding  to  given 
H,p  and  the  function  o(z)  on  tO,H)? 

According  to  the  basic  notion  recalled  at  the  end  of 
SI,  a  progressively  larger  flow  force  S  must  be  manifested 
in  the  channel  as  the  process  of  extraction  causing  the  flow 
is  intensified.  But  the  results  of  §52  and  3  show  that  the 


maximum  possible  S  achievable  without  the  entrainment  of 

fluid  originally  above  z  -  H  is  realized  by  the  critical  self- 

2 

similar  flow  with  y  -  jz.  This  flow  therefore  comprises  an 
upper  limit  for  the  possible  drawdown  condition.  If  drawdown 
does  not  occur  until  increased  extraction  raises  S  to  S  , 

Itl 

then  this  flow  is  necessarily  realised  at  the  limit,  and  no 
further  increase  of  S  is  possible  without  drawdown. 

2.  The  principle  of  maximum  flow  force  was  demonstrated 
in  SS2  and  3  without  regard  to  the  grading  of  the  contraction 
through  which  the  flow  approaches  the  channel,  but  this  factor 
evidently  will  determine  whether  the  limiting  flow  is  realizable. 
It  is  plausible  that  this  flow  does  precede  drawdown  when 
the  contraction  is  extremely  gradual.  Otherwise,  as  the  known 
behaviour  of  open-channel  flows  suggests,  an  approach  to  the 
critical  condition  is  liable  to  be  interrupted  by  wave 
formation.  [Note  that  in  practice  it  is  found  difficult  to 
produce  a  smooth  open-channel  flow  in  a  slightly  subcritical 

lc 

condition,  say  with  0.6  <  F  <  1,  where  F  =  u/(gh) 

(cf.  Binnie  et  al.  1955).] 

A  wavy  flow  in  the  channel  itself  is  a  possible  precursor 
of  drawdown  when  the  contraction  is  no*-  gradual  enough,  but 
other  possibilities  are  indicated  by  the  findings  of  S5. 

Suppose  that  the  contraction  is  very  gradual  but  the  sink  of 
fluid  at  the  end  of  the  channel  is  so  arranged  that  development 
of  a  self-similar  flow  is  hindered  (e.g.  there  are  two  slots 
extracting  fluid) .  Then  as  S  is  raised  towards  the  value  at 


which  drawdown  begins,  the  flow  may  still  differ  appreciably 
from  the  self-similar  flow,  and  consequently  a  local  crisis 
giving  rise  to  waves  may  occur  some  way  upstream. 

3.  It  must  be  acknowledged  that  the  present  estimate 
of  the  drawdown  condition  disagrees  radically  with  the  view 
of  the  matter  proposed  by  Huber  (1960),  who  calculated  a 
critical  condition  at  which  the  flow  of  a  stratum  of  homogeneous 
perfect  fluid  towards  a  line  sink  first  entrains  a  superposed 
layer  of  fluid  with  smaller  density.  The  present  interpretation 
also  conflicts  in  principle,  but  is  in  its  outcome  more 
easily  reconciled,  with  a  calculation  by  Craya  (1949)  on  a 
quite  different,  approximate  basis.  The  two  estimates  were 
discussed  by  Yih  (1965,  p.128),  who  cited  unpublished 
experimental  results  showing  a  large  discrepancy  with  Huber's 
theoretical  prediction.  In  terms  of  the  Froude  number  F  for 
the  flow  in  the  lower  stratum  (with  the  density  difference 
incorporated  into  F  in  the  usual  way) ,  the  drawdown  condition 
was  calculated  to  be  F  =  1.66,  whereas  according  to  present 
ideas  it  is  just  F  =  1.  The  former  value  can  at  once  be 
rejected  as  a  practical  threshold  for  flows  originating  from 
a  large  reservoir,  because  such  flows  with  F  >  1  and  stagnant 
fluid  above  are  impossible  in  the  absence  of  a  throat  upstream  - 
a  feature  not  recognized  in  Huber's  model.  The  experimental 
drawdown  condition  reported  by  Yih  is  roughly  F  ■  0.7,  which 


is  quits  consistent  with  the  present  interpretation  when 
allowance  is  made  as  suggested  for  the  possibility  of  wave 
formation  forestalling  the  idealized  critical  condition  P  -  1. 

4.  Cases  where  p  is  continuous  at  z  ■  H  provide  a 

different  interpretation.  Note  that  whatever  the  flow  force 

determined  by  the  process  of  extraction  at  the  end  of  the 

channel,  there  is  a  least  value  of  H  for  which  the  required  S 

can  be  realized.  Then  S  ■  S  (H) ,  where  S  is  the  maximum 

m  m 

of  S  for  a  given  H,  as  given  by  (2.8).  In  other  words,  for 

a  given  flow  force,  either  a  flow  is  developed  having 
2 

y  »  jz  in  [0,H]  for  the  respective  minimum  H,  or  a  deeper 
layer  of  fluid  is  drawn  into  motion.  At  least  in  a  straight 
channel  following  a  gradual  contraction,  it  may  therefore  be 
expected  that  self-similar  flows  will  tend  to  be  realized 
whenever  fluid  is  withdrawn  from  a  continuously  stratified 
reservoir. 

5.  As  is  known  for  open-channel  flows  of  a  homogeneous 
fluid,  steady  stratif ied-fluid  flows  along  a  convergent- 
divergent  channel  are  likely  to  be  in  better  accord  with 
rfiallow-water  theory  than  near-critical  flows  in  a  straight 
channel.  An  everywhere  subcritical  flow  with  ay/fcx  »  0  at 
the  minimum  section  is  a  theoretical  possibility;  but,  as 
the  open-channel  analogy  shows,  it  is  liable  to  be  swept  away 
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in  the  divergent  part  of  the  channel  unless  the  final  outflow 
is  specially  restricted.  The  raising  of  downstream  flow  force 
concomitantly  with  increasing  the  extraction  rate  will 
generally  produce  a  flow  that  is  critical  at  the  throat  and 
supercritical  downstream.  The  principle  demonstrated  in  §6 
accordingly  indicates  that  a  self-similar  flow  will  tend  to 
be  established  in  this  situation,  since  it  gives  rise,  at 
the  downstream  end  of  the  channel,  to  a  flow  force  that  is 
the  maximum  possible  without  additional  fluid  being  extracted. 
This  conclusion  must  be  regarded  with  caution,  however,  in 
view  of  the  artificial  feature  that  a  layer  of  stagnant  fluid 
deeper  than  |h  lies  above  the  fast,  supercritical  flow.  The 
possibility  that  instabilities  of  the  Kelvin-Helmholtz  type 
may  precipitate  a  dissipative  transition  (hydraulic  jump) 
back  to  subcritical  flow,  also  the  possibility  of  a  super¬ 
critical  drawdown  condition  such  as  found  by  Huber  (1960), 
suggests  that  the  flows  in  question  may  not  be  realizable  in 
a  far  supercritical  condition. 

It  is  noteworthy  that  Wood  (1968),  considering  a  two- 
layer  model,  derived  a  self-similar  flow  in  a  convergent- 
divergent  channel  on  the  basis  of  the  hypothesis  that  the 
shallow-water  equations  should  have  a  smooth  solution  everywhere. 
His  result  agrees,  of  course,  with  the  observations  made  in 
SS4  and  5,  but  the  present  flow-force  principle,  rather  than 
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an  arbitrary  hypothesis  of  smoothness,  is  much  more  telling 
as  a  reason  why  the  self-similar  flow  should  be  generated. 

Wood  also  presented  some  experimental  results  approximately 
confirming  his  prediction. 

6.  The  present  investigation  has  focussed  on  the 
case  of  flow  in  layers  lying  on  a  horizontal  plane,  but  all 
aspects  of  the  theory  extend  more  or  less  immediately  to  the 
case  of  withdrawal  from  internal  layers.  The  first  boundary 
condition  in  (4.4)  has  to  be  replaced  by  another,  akin  to 
the  second  of  (4.4),  applying  at  the  lower  interface  with 
stagnant  fluid,  and  the  self-similar  solution  corresponding 
to  (4.5)  -  (4.7),  but  now  defined  on  [-H.H] ,  say,  is 
appropriately  modified.  The  example  considered  in  §2,  for 
instance,  extends  precisely  to  a  solution  like  (5.1)  on 
[-1,1],  and  the  even  as  well  as  odd  Legendre  polynomials  are 
then  required  to  express  an  arbitrary  perturbation  in  the 
manner  of  (5.6)  and  (5.7). 

This  work  was  completed  during  an  extended  visit  to 
the  Mathematics  Research  Centej  at  the  University  of  Wisconsin 
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colleagues  there,  particularly  Professor  John  Nohel,  Director 
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